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ABSTRACT 



Context. In quiescence, the proposed optical counterpart to the neutron star low mass X-ray binary 4U 2129+47 (V1727 Cyg) shows 
a spectrum consistent with a late F-type subgiant and no radial velocity variations on the 5.24 hour binary period. This could imply 
that V1727 Cyg is a chance line of sight interloper. Radial velocity measurements, however, showed evidence for a longer term ~40 
km s"' shift, which suggested that 4U 2129+47 could be a hierarchical triple system, with the F-type star in a wide orbit about the 
inner low mass X-ray binary. 

Aims. In order to confirm the long-term radial velocity shift reported in Garcia et al. (1989) and its amplitude, we obtained spectro- 
scopic observations of V1727 Cyg during 1996 and 1998 with the William Herschel Telescope using the ISIS spectrograph. 
Methods. We determined radial velocities from the ISIS spectra by means of the cross-correlation technique with a template spectrum. 
Results. The resulting radial velocities show variations with a maximum amplitude of ~ 40 km s"', confirming previous results and 
supporting the F-type star as being the third body in a hierarchical triple system. The odds that this star could be an interloper are 
~ 3 X 10-". 
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1. Introduction 

4U 2129+47 was discovered as an active, but weak X-ray source 
in the fourth UHURU survey (Forman et al. 1978), and its op- 
tical counterpart was identified as VI 727 Cyg (Thorstensen et 
al. 1979). V1727 Cyg was found to be of ~17th magnitude, and 
.5^ ' exhibited a large amplitude (AB ~1.5) photometric periodicity 
K^ of 5.24 hours. This was taken to be indicative of the orbital pe- 
^ riod, with the photometric modulation ascribed to X-ray heat- 
ed ing of the companion (Thorstensen et al. 1979). Analysis of the 
X-ray light curve led to the development of a model where the 
system is viewed edge on and surrounded by an accretion disc 
corona (White & Holt 1982). Attempts to measure the mass of 
the compact object produced confusing results, with radial ve- 
locity studies implying a compact object mass of 0.6 + 0.2 M© 
(Home et al. 1986). This led to the conclusion that the binary 
system was a cataclysmic variable. However, the discovery of 
an X-ray burst confirmed the identity as a neutron star X-ray bi- 
nary (Garcia & Grindlay 1987). Assuming a neutron star primary 
and a 5.24 hour period, the empirical relationship described by 
Robinson (1976) between the companion mass and the binary 
period indicates a mass of 0.59 < M(Mq) < 0.80 and a radius of 
0.59 < R{Re) < 0.68 (Thorstensen et al. 1979). These parame- 
ters suggest a K - M V spectral type for the companion star. 

X-ray observations in 1983 failed to detect 4U 2129+47, and 
contemporaneous optical observations showed that VI 727 Cyg 
had dimmed to ~18.5th mag (Pietsch et al. 1986). This period 
of quiescence was a good opportunity for detecting the compan- 
ion star with a view to constraining the orbital parameters of 



the system. The first such study was conducted by Garcia et al. 
(1989). Spectra were taken during four observing runs between 
June 1987 and October 1988 in order to obtain more accurate 
measurements of the neutron star mass without the strong influ- 
ence of accretion disc emission. The results were surprising: no 
variations on a 5.24 hour period were detected. From the upper 
limits on the orbital radial velocity variations, the neutron star's 
mass appeared to be < 0.1 M© (Garcia et al. 1989), which is ob- 
viously inconsistent with a neutron star primary. Moreover, the 
June 1987 run provided a mean velocity significantly higher (by 
~ 40 km s"') than observed in the other nights. Another sur- 
prise was the identification of V1727 Cyg as an F7-IV star - a 
stable 5.24 hour orbit about a neutron star is smaller than the ra- 
dius of an F7-IV star (by a factor of ~ 1.5). These results have 
led to the suggestion that 4U 2129+47 is a hierarchical triple 
system. The observed F-type star is postulated to be in a wide 
orbit about the centre of mass of the close pair, while the inner 
companion is perhaps a K-type dwarf (Garcia et al. 1989). The 
F7 star dominates the optical spectrum in the cuiTent quiescent 
state, and so radial velocity shifts on a 5.24 hour period are not 
observed. Garcia et al. (1989) suggested that a 30 day outer pe- 
riod could account for the on-ofF cycles seen in 4U 2129+47. 
An outer period on this order would drive a periodic eccentric- 
ity in the inner binary on a 45 year timescale, consistent with 
the on and off states that have been observed since the 1930s. 
The radial velocity data previously obtained and newly reported 
herein are consistent with such a period. Chevalier et al. (1989) 
have pointed out that even without the observed radial velocity 
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shift, the lack of large radial velocity variations would lead to 
the postulation of a third stellar component. 

An obvious alternative to the triple hypothesis is a chance 
line of sight alignment. Ground based studies show that the 
on and off state optical counterparts are coincident to 0.26" 
indicating that the likelihood of a chance alignment is 10"^ 
(Thorstensen et al. 1988). HST imaging shows that there is 
no nearby companion of comparable magnitude within 0.04" 
(Deutsch et al. 1996). Scaling this radius to the ground based re- 
sults further reduces the chance alignment probability to 2x10"^. 

In this paper, we present spectra of V 1 727 Cyg taken in order 
to analyse its systemic radial velocity over a time interval similar 
to the long (~ 30 day) period predicted for the late F-type star 
if this is the outer component of a triple system (Garcia et al. 
1989). Previous radial velocity data were acquired during short 
observing runs (over at most 2 nights), separated by periods of 
many months. In order to test the proposed ~ 30 day period it 
was necessary to conduct a study over a time interval of ~ 1 
month, with data points separated by about a week. 



2. Observations 

Optical spectra were acquired in service mode during 1996 and 
1998 using the 4.2m William Herschel Telescope on La Palma 
equipped with the dual-beam ISIS spectrograph. V 1 727 Cyg was 
observed with both arms of the instrument. 

The 1996 data set was acquired on August 5/12/19/25 UT. 
The observations were made using the R1200R grating and 
TEK2 CCD in the red arm and the R600B and TEKl CCD 
in the blue arm, giving a dispersion of 0.40 A pixel"' and 
0.78 A pixel"', respectively. The slit width used ranged from 
1" .0 to 1".3. From measurements of the arc-lamp and night 
sky emission lines, we find that this yielded a spectral resolu- 
tion from 1 .0 to 1 .2 A (red arm) and 1 .4 to 2.0 A (blue arm). The 
red and blue data cover the wavelength range AA6370 - 6760 
and /1/14340 - 4910. On each night, two to three 30 min spec- 
tra of VI 727 Cyg were taken, along with spectra of a template 
star, either BD+47 4219 or HD 222368 (both F7-IV stars). HD 
222368 is a GCRV star with a velocity of 5 .0 + 0.9 km s" '. In to- 
tal, 1 1 spectra of VI 727 Cyg were obtained. Calibration frames 
were also acquired, in particular arc lamp spectra. However, the 
night of August 4 had no flat fields taken. The data have been 
included in the analysis as the signal to noise was high enough 
to mask the fixed pattern noise: performing a cross correlation 
between an August 5 object spectra and an identically extracted 
strip from a flat field produced no peak, indicating that the fixed 
pattern noise is small. The night of August 12 (with 2 target spec- 
tra) was discarded, as the arc lamp spectra were not rigorously 
taken close in time to the target spectrum, leading to an uncertain 
wavelength calibration. 

The 1998 data set consists of 6 nights of data: June 19/27, 
July 3/17/22 and August 2 UT. The observations were taken 
with the R1200R grating and TEK2 CCD on the red arm and the 
R1200B grating and EEVIO CCD on the blue arm of the spec- 
trograph. The slit width used ranged from 1".0 to 1".2. In the red 
arm, this provided a spectral resolution from 0.8 to 1.5 A. The 
blue arm data was of poor quality and were not used in the analy- 
sis. The red data covered the wavelength interval /1/163 10 -67 10. 
During each night two to three spectra of VI 727 Cyg with inte- 
gration times ranging from 15 to 30 min were taken, along with 
all associated calibration frames. Spectra of BD-i-47 4219 were 
also acquired each night (except on July 3). A total of 15 spectra 
of VI 727 Cyg were taken. The two spectra acquired on June 27 



were discarded due to their low quality. Furthermore, the night 
of July 16 was unsuitable for use due to instrumental problems. 
Both data sets were reduced with standard IRAF routines and 
the spectra were extracted with the IRAF KPNOSLIT package. 



3. Analysis 

Velocities were computed using the Fourier cross-correlation 
technique developed by Tonry & Davis (1979) and implemented 
in the IRAF task FXCOR. A sum of three HD 222368 spectra 
(taken on 1996 August 18 when conditions were very good) was 
used as a template spectrum. Fig. 1 shows this composite spec- 
trum from the red arm along with a sum of four VI 727 Cyg spec- 
tra. Prior to the cross-correlation, the target and template spectra 
were re-sampled into a common logarithmic wavelength scale 
and normalized by dividing with the result of fitting a low-order 
spline to the continuum. Correlation was performed without Ha 
in the red, to ensure that any possible residual emission from 
the accretion disc did not interfere with the measured velocity 
of V1727 Cyg. In the blue, H/3 was included, as (besides the 
G-band) the spectra lack well defined features in the available 
region, and without H/3 correlation fits were poor. 

Fig. 2 and Table 1 show the radial velocity of V1727 Cyg 
from the 1996 and 1998 data sets. The radial velocities shown in 
Fig. 2 are weighted averages of radial velocities obtained from 
the red and (when available) blue spectra. All velocities are he- 
liocentric. Table 1 also shows the average radial velocities ob- 
tained as derived using only the blue or ned spectra. As in Garcia 
et al. (1989) the errors on the individual velocities produced by 
FXCOR have been computed as cr = C/(l + r), where C is 
a constant determined from the observed scatter in the veloc- 
ity of comparison stars, and r is the correlation coefficient (see 
Tonry and Davis 1979). Gaussian statistics were used through- 
out when calculating mean values. Checking the wavelength of 
night sky emission lines served as a zero point calibration for 
the data: these lines are seen only at wavelengths longer than 
~ 5000 A, so this method could not be used for the blue frames. 
The average night sky line velocity is approximately -0.4 km s"' , 
which is consistent with zero given the errors. The wavelength 
of each night sky line was found by finding the centroid of a 
Gaussian function fitted to each line. Both the larger error bar 
for the VI 727 Cyg radial velocity and the lack of sky line mea- 
surements on 1998 July 3 are due to the fact that the spectra were 
obtained during dawn. As noted in section 2, no radial velocity 
templates were observed during this night. 



4. Discussion 

Our data do not show the expected short term, ~ 300 km s"' 
semi-amplitude radial velocity variations indicative of a neu- 
tron star primary. This is consistent with all previous stud- 
ies of V1727 Cyg in quiescence (see e.g. Garcia et al. 1989; 
Garcia 1989; ChevaUer et al. 1989; Cowley & Schmidtke 1990). 
Furthermore, the longer term velocity variations found by Garcia 
et al. 1989 are confirmed in this new dataset. The 40 km s"' am- 



plitude of the variation is similar to that found earlier, and the 
variation is detected between 7 and 10 cr level. 

If we accept that the late F-type subgiant is indeed orbiting 
an unseen companion (or companions) we can ask what are the 
odds that it could still be a line of sight interloper, rather than 
physically associated with 4U 2129+47? The fraction of multi- 
ple systems in the Galaxy is an active topic of study, with esti- 
mates ranging from 10% to 80% (Lada 2006). It is now clear that 
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the fraction depends upon mass, with more massive stars more 
likely to be members of multiple systems. Also, the majority of 
the systems are widely detached and show only small (few km 
s~') radial velocity variations. We used the study of Duquennoy 
& Mayor (1991) to estimate the fraction of F7 stars that would 
show velocity variations of at least 40 km s"' (corresponding to a 
radial velocity semi-amplitude K of 20 km s"'). This study mea- 
sured orbital velocity variations of a sample of F7 to G8 stars 
in order to determine orbital periods and multiplicity frequency. 
The mean period in the sample was logiP) - 4.8 with a roughly 
Gaussian dispersion of log(P) = 2.3 with P in days. From Table 
2 of Duquennoy & Mayor (1991) we estimate that K = 20 km 
s"' corresponds approximately to P = 100 days. The fraction of 
systems with periods less than 100 days, and therefore likely to 
show velocity variation of 40 km s"' or more, is approximately 
1 in 7. Therefore the odds that VI 727 Cyg could be an interloper 
are reduced from 2 x lO"*" to approximately 3 x lO"*". 

Recently, Bozzo et al. (2007) have presented time analysis 
of two XMM-Newton observations of 4U 2129-1-47 in quies- 
cence obtained ~ 22 days apart, finding evidence for a delay 
of ~ 200 s for the mid-eclipse times measured from the X-ray 
observations. This delay can be explained as due to the orbital 
motion of the compact 4U 2129-1-47 binary around the center 
of mass of a triple system. In light of these new observations, 
coupled with the radial velocity measurements presented in this 
paper, we conclude that the triple hypothesis best explains all 
features of 4U 2129-1-47, and that the late F-type subgiant is the 
outer component of the triple system. 



5. Conclusions 

The X-ray source 4U 2129-1-47 has been a candidate hierarchi- 
cal triple for some time. The optical counterpart, V1727 Cyg, 
has features inconsistent with those expected for a close binary 
companion, and previously displayed radial velocity variation 
over several weeks. The new spectroscopic measurements pre- 
sented here confirm this long term shift of ~ 40 km s ', which 
strongly suggests that 4U 2129-1-47 is indeed a triple. However, 
the outer period is yet to be observed directly. Doing so would 
allow for more sophisticated modeling of the system. It may also 
help understand the evolutionary history of 4U 2129H-47, which 
as of now is somewhat elusive. 
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Table 1. Heliocentric radial velocities in km s ') for V1727 Cyg 
in quiescence. 



Date Radial Velocity Radial Velocity Radial Velocity 

(UT) (blue and Red arm) (blue arm) (red arm) 



1996 Aug 5.0 


-48.4 ±5.8 


-48.7 ±7.8 


-48.3 ± 2.7 


1996 Aug 19.2 


-45.4 ±5.1 


-48.3 ± 6.0 


-39.3 ± 34 


1996 Aug 25.1 


-8.0 ±3.9 


-3.1 ±4.2 


-144 ±4.0 


1998 Jun 19.2 


-40.3 ± 3.6 


— 


-40.3 ± 3.6 


1998 Jul 3.2 


-5 ±20 


— 


-5 ±20 


1998 Jul 22.1 


-17.2 ±5.2 


— 


-17.2 ±5.2 


1998 Aug 2.1 


-22 ±11 


— 


-22 ±11 
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Fig. 1. Summed spectra of Willi Cyg compared to that of the 
template spectrum used, HD222368. 
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